Closed orbits and their bifurcations in the crossed-fields hydrogen atom 
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A systematic study of closed classical orbits of the hydrogen atom in crossed electric and magnetic 
fields is presented. We develop a local bifurcation theory for closed orbits which is analogous to 
the well-known bifurcation theory for periodic orbits and allows identifying the generic closed-orbit 
bifurcations of codimension one. Several bifurcation scenarios are described in detail. They are 
shown to have as their constituents the generic codimension-one bifurcations, which combine into 
a rich variety of complicated scenarios. We propose heuristic criteria for a classification of closed 
orbits that can serve to systematize the complex set of orbits. 
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I. INTRODUCTION 

Closed-orbit theory |TJ ^ has proven to be the key tool 
to analyze the photoabsorption spectra of atoms in ex- 
ternal fields, ft interprets spectral oscillations semiclassi- 
cally in terms of closed orbits of the underlying classical 
system, i.e. of classical orbits starting at and return- 
ing to the nucleus. A complete semiclassical description 
of an atomic spectrum therefore requires a sufficiently 
detailed understanding of the classical closed orbits. In 
particular, the possible types of closed-orbit bifurcations 
must be described, so that the generation of new closed 
orbits upon varying the external field strengths can be 
accounted for. 

For the hydrogen atom in a magnetic field, the sys- 
tematics of closed orbits and their bifurcations has been 
known for a long time |, |, |, |, @, §. For the hy- 
drogen atom in crossed electric and magnetic fields, the 
classical mechanics is much more complicated because 
three non-separable degrees of freedom have to be dealt 
with. Although a number of closed orbits have been 
identified in experimental or theoretical quantum spec- 
tra (|, ||, |l0|, p] |lj] , a systematic study of these orbits 
and their bifurcations is still lacking. 

Considerable effort has been spent during the past 
decade on the study of the classical mechanics of the 
crossed-fields hydrogen atom in the limit of weak external 

fields 0. |53 HI. 0- Tliat work relies 011 the observa- 
tion that for weak external fields the principal quantum 
number n, or its classical analogue n = is con- 

served to a higher degree of precision than the angular- 
momentum and Lenz vectors L and A. The latter are 
conserved in the pure Kepler problem, but acquire a slow 
time-dependence in weak fields, so that the electron can 
be visualized as moving on a slowly precessing Kepler 
ellipse. 

The most important result in the present context is the 
finding first described in [Q that there are four Kepler 
ellipses that remain unperturbed by the external fields 
to first order in the field strength, i.e. among the con- 
tinuous infinity of periodic orbits of the unperturbed Ke- 
pler proble there are four orbits that are periodic even in 
the presence of external fields. These fundamental peri- 



odic orbits can be regarded as the roots of "family trees" 
of periodic orbits. More complicated orbits are created 
out of the fundamental orbits by bifurcations as the field 
strengths increase. 

However, none of the fundamental periodic orbits is 
closed at the nucleus. Their knowledge therefore does 
not aid in the classification of closed orbits. A first sys- 
tematic study of closed orbits in the crossed-fields system 
and their bifurcations was performed by Wang and Delos 
fL8| . These authors presented orderly sequences of bifur- 
cations of planar closed orbits (i.e. orbits in the plane 
perpendicular to the magnetic field), which they inter- 
preted in terms of an integrable model Hamiltonian. 

In the present paper we undertake a systematic in- 
vestigation of closed orbits and their bifurcations in the 
crossed-fields hydrogen atom. In section B the sym- 
metries of the Hamiltonian are briefly reviewed. Sec- 
tion HI A presents the general framework of a lo cal bifur- 



cation theory of closed orbits, and section III B describes 
the generic codimension-one bifurcations. A discussion 
of complex ghost orbits is included in each case because 
they are known to play an important role in semiclassics 
p9j |20f . In section IV, the closed orbits in the hydro- 
gen atom in a magnetic field are surveyed. Section ^ 
then details the bifurcation scenarios actually observed 
in the crossed-fields system. It is shown that the ele- 
mentary codimension-one bifurcations actually form the 
building blocks of the bifurcations scenarios, but that in 
many cases complicated scenarios consisting of several el- 



ementary bifurcations occur. In section VI, a heuristic 



classification scheme for the closed orbits in crossed fields 
is proposed which is based on the well-known classifica- 
tion for the closed orbits in a magnetic field. The actual 
semiclassical quantization of the crossed-fields hydrogen 
atom in the framework of closed-orbit theory, which is 
based on the results presented here, is described in an 
accompanying paper [Ell . 



II. THE CLASSICAL HAMILTONIAN 

Throughout this work, we will assume the magnetic 
field to be directed along the z-axis and the electric field 
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TABLE I: The symmetry transformations of the crossed-fields 
system: Transformation of initial and final angles and sym- 
metry conditions for doublets. Singlets satisfy #j = t?f = ^ 
and cpi = — iff. 

to be directed along the x-axis. In atomic units, the 
Hamiltonian describing the motion of the atomic electron 
then reads 

H = \p 2 - \ + \BL Z + X -B 2 p 2 + Fx , (1) 

where r 2 = x 2 + y 2 + z 2 , p 2 = x 2 + y 2 and L z is the z- 
component of the angular momentum vector. By virtue 
of the scaling properties of the Hamiltonian ([!]), the dy- 
namics does not depend on the energy E and the field 
strengths B and F separately, but only on the scaled en- 
ergy E = B~ 2 / Z E and the scaled electric field strength 
F = B~ 4 ^ 3 F. Upon scaling, all classical quantities 
are multiplied by suitable powers of the magnetic field 
strength B. In particular, classical actions scale accord- 
ing to S = B _1 / 3 5. These scaling prescriptions will be 
used throughout this work. 

In crossed fields, two angles are required to characterize 
the starting or returning direction of a closed orbit. We 
will use the polar angle i? between the trajectory and the 
magnetic field axis and the azimuthal angle <p between 
the projection of the trajectory into the x-y-plane and 
the electric field axis. 

The hydrogen atom in crossed fields does not possess 
any continuous symmetries so that, apart from the en- 
ergy, no constant of the motion exists and three non- 
separable degrees of freedom have to be dealt with. There 
are, however, three discrete symmetry transformations of 
the crossed-fields system, namely 

• the reflection Z at the x-y-plane, 

• the combination T of time-reversal and a reflection 
at the x-z-plane, 

• the combination C=ZT of the above. 

The effects of the transformations on the initial and final 
angles of a closed orbit are summarized in table [I]. 

The application of these transformations to a given 
closed orbit yields a group of four closed orbits of equal 
length. Typically, these orbits will all be distinct, so that 
closed orbits in the crossed-fields system occur in quar- 
tets. In particular closed orbit can be invari- 
ant under one of the symmetry transformations. In this 
case, there are only two distinct orbits related by symme- 
try transformations. We will refer to them as a doublet, 



or more specifically as a Z-, T-, or C-doublet, giving the 
transformation under which the orbits are invariant. The 
transformations of the initial and final angles listed in ta- 
ble I yield symmetry conditions that an orbit invariant 
under any of the transformations must satisfy. These are 
also given in table |. In special cases, a closed orbit can 
be invariant under all three symmetry transformations. 
It then occurs as a singlet, since no distinct orbits can be 
generated from it by a symmetry transformation. 

Among the symmetry transformations, the reflection Z 
plays a special role in that it is a purely geometric trans- 
formation. There is, therefore, an invariant subspace of 
the phase space, viz. the x-y-plane perpendicular to the 
magnetic field. This plane is invariant under the dynam- 
ics and therefore forms a subsystem with two degrees of 
freedom. 

In connection with bifurcations of orbits it is essential 
for semiclassical applications to study complex "ghost" 
orbits along with the real orbits, i.e. to allow coordinates 
and momenta to assume complex values. For ghost or- 
bits, another reflection symmetry arises, viz. the sym- 
metry with respect to complex conjugation. Since the 
Hamiltonian (|l|) is real, it is invariant under complex con- 
jugation. Therefore, ghost orbits always occur in pairs of 
conjugate orbits. 



III. CLOSED-ORBIT BIFURCATION THEORY 
A. General theory 

The dynamics of the hydrogen atom in a pure magnetic 
field possesses time-reversal invariance if it is restricted to 
the subspace of vanishing angular momentum L z . Due 
to this symmetry, an electron returning to the nucleus 
will rebound from the Coulomb center into its direction 
of incidence and retrace its previous trajectory back to 
its starting direction. Therefore, any closed orbit is ei- 
ther itself periodic or it is one half of a periodic orbit. 
Due to the close link between closed orbits and periodic 
orbits, closed-orbit bifurcations can be described in the 
framework of periodic-orbit bifurcation theory developed 
by Mayer ||, £2|. In particular, in a magnetic field closed 
orbits possess repetitions, so that arbitrary m-tupling bi- 
furcations are possible. 

In the presence of crossed electric and magnetic fields, 
the time-reversal invariance is broken, and no general 
connection between closed orbits and periodic orbits re- 
mains. As a consequence, the techniques of periodic-orbit 
bifurcation theory are no longer applicable, and a novel 
approach to the classification of closed-orbit bifurcations 
must be found. In this section, a general framework for 
the discussion of closed-orbit bifurcations will be intro- 
duced. 

The crucial step in the development of the bifurcation 
theory of periodic orbits is the introduction of a Poincare 
surface of section map in the neighborhood of the orbit. 
The Poincare map describes the dynamics of the degrees 
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of freedom transverse to the orbit, and the orbit bifur- 
cates when the transverse dynamics becomes resonant 
with the motion along the orbit. 

For periodic orbits, a Poincare map is specified by fix- 
ing a surface of section in phase space which is transverse 
to the orbit. For a point P on the surface of section, the 
trajectory starting at P is followed until it intersects the 
surface of section again. This intersection point is de- 
fined to be the image of P under the Poincare map. The 
periodic orbit itself returns to its starting point, so that 
it appears as a fixed point of the Poincare map. 

This prescription is not directly applicable to closed 
orbits because they do not return to their starting point 
in phase space. Therefore, a trajectory starting on the 
surface of section will not in general intersect the surface 
again. To arrive at a meaningful definition of a Poincare 
map, one must use two surfaces of section: the first trans- 
verse to the initial direction of the orbit, the second trans- 
verse to its final direction. A trajectory starting in the 
neighborhood of the closed orbit on the initial surface 
of section £j will then have an intersection with the fi- 
nal section £/, so that a Poincare map is well defined. 
As in the case of a periodic orbit, the Poincare map is 
symplectic. 

Unlike with periodic orbits, the notion of a closed or- 
bit is not invariant under canonical transformations. The 
distinction between position space and momentum space 
must therefore be kept. Let {qi,p%) and (q/,Pf) be canon- 
ical coordinates on the surfaces £j and E f chosen so that 
qi and qf are position space coordinates in the directions 
perpendicular to the initial or final directions of the orbit. 
The origins of the coordinate systems are fixed so that the 
position of the nucleus is qi = or qf — 0, respectively. 
Closed orbits are then characterized by qi — qf = 0. 
In crossed fields three spatial dimensions must be dealt 
with, so that each of qi,Pi,qf,Pf is a two-dimensional 
vector. The reader may conveniently picture qi and qf 
as Cartesian coordinates, although in this case the con- 
jugate momenta pi and pf diverge as the Coulomb singu- 
larity is approached. This difficulty can be overcome by 
means of a Kustaanheimo-Stiefel regularization p3[ . Co- 
ordinates having the properties described above can then 
be shown to exist, as will be discussed in detail elsewhere 

A closed orbit can start in £j with arbitrary initial 
momentum p.;, but it must start in the plane qi — 0. The 
Poincare map maps this plane into a Lagrangian manifold 
in £/. Closed orbits are given by the intersections of this 
manifold with the plane qf = 0. In a less geometrical way 
of speaking, closed orbits can be described as solutions of 
the equation qf(pi,qi — 0) = 0. A particular solution of 
this equation, corresponding to the orbit the construction 
started with, is given by qf(j>% — 0) = 0. If the matrix 
B = dqf/dpi is non-singular at pi = 0, this solution 
is locally unique and, by the implicit function theorem, 
will persist upon the variation of parameters. Thus, the 
closed orbit cannot undergo a bifurcation unless M = 
det B = 0. 



type 


transformation 


regular matrix 


Fi(qi,qf) 


Pi 


= +dF 1 /dg t , P f = -dFx/dqj 


B 


F 2 (qi,Pf) 


Pi 


= +dF 2 /dqi , qf = +dF 2 /dpf 


D 


F 3 (pi,qf) 


qi 


= -dFs/dpi , p f = -dF 3 /dq f 


A 




qi 


= -dFi/dpi , q f = +dF 4 /d Pf 


C 



TABLE II: Overview of generating functions of different types 
(cf. @). 



An overview of the bifurcation scenarios to be expected 
when det B — can be obtained from a description of the 
possible modes of behavior of the Poincare map. This 
can most conveniently be achieved if the Poincare map is 
represented by a generating function [ prj . The generating 
function can be chosen to depend on any combination 
of initial and final positions and momenta, as long as 
they form a complete set of independent coordinates. We 
adopt the well-known conventions of Goldstein |2(| for 
denoting different types of generating functions, which 
are summarized in table 0. 

For a generic symplectic map, all possible sets of coor- 
dinates and momenta are independent, so that generating 
functions of any type exist. At a closed-orbit bifurcation, 
however, a degeneracy indicated by the condition that 
B = dqf/dpi be singular arises, so that care must be 
taken in choosing a generating function. Loosely speak- 
ing, if B is singular, pi cannot be determined from 
and qt, so that it may be conjectured that no generating 
function of type F\ exists. To confirm this conjecture, 
we study a linear symplectic map 



q f = Aq l + Bpi 



Pf =Cqi + Dpi 



(2) 



with four matrices A, B, C, D satisfying the symplecticity 
conditions |27| 

A T C = C T A , B T D= D T B , A T D - C T B = 1 , 

AB T = BA T , CD T = DC T , AD T - BC T = 1 , 

(3) 

where T denotes the transpose. A generating function 
for the linear map (|^) must be quadratic in its variables. 
From the ansatz 

Fife, q f ) = \q]Rq S + q] S Ql + \qjTq, (4) 

with matrices R, S, T, the map (Q) is obtained if 

R=-DB- 1 , T=-B~ 1T A, 
S = B- 1T = DB^A-C . 



(5) 



The two expressions given for S are equal by virtue of (^) . 
As expected, a generating function of type F\ does not 
exist if B is singular. A similar calculation can be car- 
ried out for the other types of generating functions. For 
each type, one of the matrices A, B, C, D must be non- 
singular. These results are given in tabic 0. Locally, 
they can be extended to non-linear maps by means of 
the implicit function theorem. 
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(a) 



Pf 



(b) 



Pf 



FIG. 1: Schematic plot of the Lagrangian manifold qi — in 
E/ for the case (a) B = dqj /dpi — and (b) D — dpf/dpi = 
0. The dashed lines indicate the position of the manifold at 
slightly varied parameter values. Only in case (a) additional 
intersections with the plane qf = can arise. 



Thus, at a bifurcation of closed orbits the Poincare 
map generically possesses generating functions of all 
types except Pi. The most convenient choice is a func- 
tion of type Fi(pi,pf). The transformation equations 
associated with this type of generating function read 



Qi 



dpi 



Qf 



dp f 



(6) 



Closed orbits are characterized by qi — qf = 0. They 
therefore agree with the stationary points of the P4 func- 
tion. The classification problem of closed-orbit bifurca- 
tion theory can thus be rephrased as the problem to de- 
termine how stationary points of a real function change 
upon the variation of parameters. This question is the 
subject of catastrophe theory |^8|, ^9], pOfl . 

Catastrophe theory studies smooth real-valued func- 
tions f(x) and f{x) defined in a neighborhood of the 
origin in an n-dimensional configuration space. They are 
said to be equivalent if there is a diffeomorphism ip(x) of 
the configuration space so that 



f(x) = f(Mx)) 



(7) 



The coordinate transformation tp maps the stationary 
points of / to those of /. In this sense, the distribu- 
tions of stationary points of / and / agree qualitatively. 
Without loss of generality it can be assumed that / and 
/ have stationary points at the origin, because any sta- 
tionary point can be moved there by a coordinate trans- 
formation. After adding a constant, one has /(0) = 0. 

/ is said to be structurally stable if any small pertur- 
bation / of / (i.e. f(x) — f(x) + eg(x) with a smooth 
function g(x) and sufficiently small e) is equivalent to /. 
Notice that catastrophe theory is a purely local theory. 
It is concerned with the structural stability or instability 
of a single stationary point and the pattern of stationary 
points that can be generated from a structurally unstable 
stationary point by a small perturbation. 

In the present context, non-bifurcating closed orbits 
correspond to structurally stable stationary points of P4, 



because a small variation of parameters will bring about 
a variation of P4 which is small in the above sense and 
preserves the stationary point. The most elementary re- 
sult of catastrophe theory states that a stationary point 
of a function / is structurally stable if its Hessian matrix, 
i.e. the matrix of second derivatives of /, is non-singular. 
For the linear symplectic transformation (^|), the P4 gen- 
erating function is 



1 



F^-pjAC^pf-pJC- 1 



~Pi + IpJC^Dpi , (8) 



so that its Hessian determinant at pi = pf = can be 
found to be 



dct Hess Fa 



det B det D 
detC 



(9) 



The Hessian matrix of P4 is thus singular if cither B = 
dqf jdpi or D — dp/ jdpi is. It has been shown above that 
bifurcations of closed orbits can only occur if dct B = 
0, i.e. a bifurcating orbit corresponds to a degenerate 
stationary point of F4 . The case det D = also leads 
to a degeneracy of P4, but it cannot be associated with 
a closed-orbit bifurcation. This can also be understood 
geometrically: As illustrated in figure [l], if det B — 0, the 
Lagrangian manifold given by qi — is tangent to the 
plane qf — 0, so that it can develop further intersections 
with that plane upon a small variation of parameters. If 
det-D = 0, the manifold is tangent to the plane pt = 0, 
whence, upon a variation of parameters, it can acquire 
additional intersections with that plane, but not with the 
plane qf — 0, so that no bifurcation of closed orbits can 
arise. 

The discussion of stationary points with degenerate 
Hessian matrices, also called "catastrophes", is simpli- 
fied considerably by the splitting lemma of catastrophe 
theory {kJ . It states that if the dimension of the config- 
uration space is n and a function / on the configuration 
space has a stationary point at the origin whose Hessian 
matrix has rank n — to, a coordinate system x±, . . . ,x n 
can be introduced in a neighborhood of the stationary 
point so that 



f(xi, ...,x n ) = g(xi, x m ) + q(x m+ i 



,*»), (10) 



where q is a non-degenerate quadratic form of n — m 
variables and the function g has a stationary point with 
zero Hessian matrix at the origin. As the non-degenerate 
stationary point of q is structurally stable, the behavior 
of the stationary points of / under a small perturbation 
is determined by g only. The number of relevant vari- 
ables is thus only to, which is called the corank of the 
catastrophe. It will be assumed henceforth that a split- 
ting according to ( |Io| ) has been carried out and the non- 
degenerate part q is ignored, so that the Hessian matrix 
of / vanishes at the origin. 

Under a small perturbation of the function /, a degen- 
erate stationary point will in general split into several 
distinct stationary points. This process will be used to 
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model bifurcations of closed orbits. The degenerate sta- 
tionary points relevant to bifurcation theory are those of 
finite codimension, i.e. those for which there are smooth 
functions gi(x), . . . , gk{x) so that any small perturbation 
of / is equivalent to 

F{x) = f(x) + a l9l {x) + ■■■ + a k9k {x) (11) 

with suitably chosen constants a,. The function F(x) 
is called an unfolding of f(x), because the degenerate 
stationary point of / can be regarded as a set of several 
stationary points that accidentally coincide and are "un- 
folded" by the parameters on. The smallest value of k 
that can be chosen in ([ll]) is called the codimension of /. 
An unfolding of / with k equal to the codimension of / 
is referred to as universal. 

In the bifurcation problem, the generating function F4 
depends on external control parameters p\ , . . . , pi such 
as, e.g., the energy E or the external field strengths. If, 
for a critical value of the parameters, F4 has a degenerate 
stationary point equivalent to that of /, in a neighbor- 
hood of the critical value F4 is equivalent to the unfolding 
( |IT| ) , where the unfolding parameters on are smooth func- 
tions of the control parameters pj. The critical parame- 
ter values themselves are characterized by the condition 
that all unfolding parameters vanish, i.e. by the set of 
equations 

ai(pi, ■ ..,pi) = Q , 

(12) 

a k (px, ...,pi) = . 

These are k equations in I unknowns. They can "gener- 
ically" only be expected to possess a solution if k < I, 
that is, the codimension of the degenerate stationary 
point must not be larger than the number of external pa- 
rameters. This construction introduces a notion of codi- 
mension for bifurcations of closed orbits which is entirely 
analogous to the codimension of bifurcations of periodic 
orbits: Bifurcations of a codimension higher than the 
number of external parameters cannot be expected to 
occur because they are structurally unstable. Under a 
small perturbation of the system they would split into a 
sequence of "generic" bifurcations of lower codimensions. 

B. Codimension-one generic bifurcations 

The considerations of the preceding section reduce the 
bifurcation theory for closed orbits to the problem of de- 
termining all catastrophes having a codimension smaller 
than the number of external parameters. In particular, 
it explains why only catastrophes of finite codimension 
are relevant. In the crossed-fields system, the number 
of parameters is two, if the scaling properties are taken 
into account. However, we will only describe bifurcations 
of codimension one in the following. They suffice to de- 
scribe the bifurcations encountered if a single parameter 
is varied while the second is kept fixed. They also give 



a good impression of the codimension-two scenarios be- 
cause a bifurcation of codimension two must split into a 
sequence of codimension-one bifurcations as soon as any 
of the parameters is changed. 

For generic functions without special symmetries, a 
list of catastrophes of codimensions up to six with their 
universal unfoldings is readily available in the literature 
pq, ps| , [30) . The classification of closed-orbit bifurcations 
presented here relies on these results. 

1. The tangent bifurcation 

There is a single catastrophe of codimension one, which 
has corank one and is known as the fold catastrophe. Its 
universal unfolding is given by 

$«(*) = ^ 3 - at , (13) 

with a denoting the unfolding parameter. The fold has 
two stationary points at 

t = ±Vo, (14) 
where it assumes the stationary values 

$ Q (±^)=T^a 3/2 . (15) 
The second derivative in the stationary points is 

$2(±Va) = ±2V^ • (16) 

The stationary points are real if a > 0. If a < 0, there 
are no stationary points on the real axis, because the 
solutions ( |l4| ) are imaginary. These complex stationary 
points correspond to closed ghost orbits in the complex- 
ified phase space. As a is varied, a tangent bifurcation 
occurs at a — 0, where two complex conjugate ghost or- 
bits turn into two real orbits or vice versa. 

All qualitative features of the bifurcation are described 
by the normal form dl3|). The stationary points, i.e. the 
closed orbits, initially move apart as \fa. A more de- 
tailed connection between the properties of the normal 
form and the closed orbits can be made in the context 
of uniform semiclassical approximations. The difference 
between the stationary values gives the action difference 
between the closed orbits, whereas the second derivatives 
- or, if the normal form has corank greater than one, the 
Hessian determinants - at the stationary points are pro- 
portional to a parameter M describing the stability of 
the closed orbit (see [^j], ^3 for details) . All these quan- 
tities are shown in figure |2j. When they are compared 
to the corresponding quantities calculated for an actual 
bifurcation in section the qualitative agreement will 
become clear. 

The fold catastrophe (|l^) describes the generation of 
two closed orbits in a tangent bifurcation. As this is the 
only generic catastrophe of codimension one, it follows 
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FIG. 2: The positions of stationary points, stationary values 
and second derivatives in the fold catastrophe. Solid lines in- 
dicate real stationary points, dashed lines complex stationary 
points. Dotted lines are coordinate axes. 

that the tangent bifurcation is the only possible type of 
closed-orbit bifurcations. In particular, once it has been 
generated a closed orbit cannot split into several orbits, 
as periodic orbits typically do. However, this statement 
needs some modification due to the presence of reflection 
symmetries in the crossed-fields system. 

2. The pitchfork bifurcation 

If the orbit under study is symmetric under one of the 
reflections, i.e. it is a singlet or a doublet orbit, the gener- 
ating function F4 in the neighborhood of this orbit must 
also possess this reflection symmetry. By this constraint, 
several of the elementary catastrophes are excluded alto- 
gether. For others, the codimension is reduced because 
the unfolding can only contain symmetric terms. 

One additional catastrophe of codimension one arises, 
viz. the symmetrized version of the cusp catastrophe 

<S> a (t) = ~t*-~at 2 . (17) 

This normal form possesses the reflection symmetry t 1— > 
— t, so that the origin is mapped onto itself under the 
symmetry transformation. There is a stationary point at 
the origin for all values of the parameter a, corresponding 
to a closed orbit which is invariant under the reflection. 
Additional stationary points are located at 

* = ±Va . (18) 

They are real if a > and are mapped onto each other 
under a reflection. Thus, the symmetric cusp (0) de- 
scribes a pitchfork bifurcation at a = 0, where two asym- 
metric orbits bifurcate off a symmetric orbit, generating 
a quartet from a doublet or a doublet from a singlet. 



FIG. 3: The positions of stationary points, stationary values 
and second derivatives in the cusp catastrophe. Solid lines in- 
dicate real stationary points, dashed lines complex stationary 
points. Dotted lines are coordinate axes. Note that station- 
ary values and second derivatives are real even for complex 
stationary points. 



The stationary values at the asymmetric stationary 
points are given by 

$ a (±V^) = -\a 2 , (19) 
the second derivative is 

$£(±Vo) = 2a . (20) 

Both the stationary values and the values of the second 
derivative are real even for a < 0, when the stationary 
points themselves are complex. Therefore, these station- 
ary points correspond to ghost orbits having real actions 
and stability determinants. The existence of this remark- 
able type of ghost orbits is again a consequence of the 
reflection symmetry: As the stationary points ( |l8| ) are 
imaginary, the reflection t 1— ► — t changes a stationary 
point and its stationary value into their complex con- 
jugates. On the other hand, the stationary values are 
invariant under the reflection, so they must be real. A 
ghost orbit having this symmetry property will be re- 
ferred to as a symmetric ghost orbit. 

The characteristic quantities of the symmetric cusp 
catastrophe are shown in figure | as a function of a. 
Again, they describe the qualitative behavior of the bifur- 
cating orbits close to the bifurcation. It should be noted 
that the stationary values (|l9| ) are negative for all val- 
ues of a, so that for a bifurcation described by fll7|), the 
actions of the asymmetric orbits must be smaller than 
those of the symmetric orbit. An alternative bifurcation 
scenario is described by the dual cusp, viz. the negative of 
dl7|). The dual cusp is inequivalent to the regular cusp, 
but the scenario it describes agrees with the above ex- 
cept that the stationary values and the second derivatives 
change their signs, so that the actions of the asymmetric 
orbits are now larger than that of the symmetric orbit. 



* 1 1 1 1 1 1 1 * 

0.2 0.4 0.6 0.8 1 1.2 1.4 

FIG. 4: Scaled actions 5* as functions of the starting angles 
■&i of closed orbits in the DKP for E = -1.4. 



IV. CLOSED ORBITS IN THE DIAMAGNETIC 

KEPLER PROBLEM 

As a basis for the description of closed orbits in the 
crossed-fields hydrogen atom, we will choose the closed 
orbits in the diamagnetic Kepler problem (DKP), i.e. in 
the hydrogen atom in a pure magnetic field. For these 
orbits a complete classification is available || ||, |[ [| 0, 
H . It will now be recapitulated briefly. 

For low scaled energies E — ► — oo, there are two funda- 
mental closed orbits: In one case, the electron leaves the 
nucleus parallel to the magnetic field until the Coulomb 
attraction forces it back. This orbit is purely Coulombic 
because the electron does not feel a Lorentz force when 
moving parallel to the magnetic field. The second closed 
orbit lies in the plane perpendicular to the magnetic field. 
Its shape is determined by the combined influences of the 
Coulomb and magnetic fields. 

Due to time-reversal invariance, both elementary or- 
bits possess arbitrary repetitions. As the scaled energy 
increases, each repetition of an elementary orbit under- 
goes a sequence of bifurcations labelled by an integer 
v = 1, 2, 3, . . . in order of increasing bifurcation energy. 
The orbits born in these bifurcations can be character- 
ized by the repetition number /z of the bifurcating orbit 
and the bifurcation number v. They are referred to || 
as vibrators if they bifurcate out of the orbit parallel 
to the magnetic field and as rotators if they bifurcate 
out of the orbit perpendicular to B. 

Further bifurcations create additional orbits from the 

V. f and R v n or "exotic" orbits not related to one of the 
two fundamental orbits. These orbits are of importance 
at scaled energies higher than those considered in this 
work, so that they will not be discussed further. For the 
scaled energy E = —1.4, the scaled actions and starting 
angles of the closed orbits are presented in figure |^. It 
can be seen that only orbits fitting into the classification 
scheme described above are present. Furthermore, orbits 
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FIG. 5: Scaled actions S and azimuthal starting angles ipi for 
planar orbits at E = —1.4 and F = 0.03 (+ symbols) and 
F — 0.05 (x symbols). 

having a common bifurcation number v lie on a smooth 
curve in the plot. For this reason, we will refer to or- 
bits characterized by a fixed v as a series of rotators or 
vibrators, respectively, and call v the series number. 



V. CLOSED ORBIT BIFURCATION 
SCENARIOS 

In the presence of a pure magnetic field, the atomic 
system possesses a rotational symmetry around the field 
axis. As a consequence, all closed orbits except for the 
orbit parallel to the magnetic field occur in continuous 
one-parameter families obtained by rotating a single or- 
bit around the symmetry axis. When a perpendicular 
electric field is added, the rotational symmetry is bro- 
ken. Out of each family, only two orbits survive [ pl| , 
or, in other words, each family of orbits splits into two 
independent orbits. 



A. Planar orbits 

The splitting of a family of orbits upon the introduc- 
tion of an electric field can most clearly be seen for planar 
orbits, i.e. for orbits lying in the plane perpendicular to 
the magnetic field. Due to the Z-symmetry, this plane is 
invariant under the dynamics. Thus, the initial direction 
of an orbit can be specified by means of the azimuthal 
angle ifi only. 

Figure || shows the actions and initial directions of the 
planar orbits for a scaled energy of E — —1.4 and scaled 
electric field strengths F = 0.03 and F = 0.05. At F = 0, 
the orbits bifurcate off a certain repetition of the planar 
closed orbit of the diamagnetic Kepler problem. For low 
F they can therefore be assigned a repetition number. It 
can clearly be discerned in figure @ from the actions of 



8 





- 


- 




(1-1 


- 
















m-4 


VJ 


















H=8 




H=10 





-0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.E 



FIG. 6: Elementary planar closed orbits at E = —1.4 and F = 
0.05. n is the repetition number. The orbits are symmetric 
with respect to the z-axis, in particular, tpi = — ipf. 



the orbits. 

As expected from the theory of the rotational symme- 
try breaking |3l| , there are two orbits for each repetition 
number, and they start in opposite directions from the 
nucleus. Moreover, the starting angle varies linearly with 
the repetition number. These findings are illustrated in 
figure |g, where for a few low repetition numbers one of 
the two orbits is shown. It can be seen that the orbits 
consist of more and more "loops" and that the starting 
angle increases regularly. The shapes are symmetric with 
respect to the x-axis because the orbits are invariant un- 
der the T-transformation, i.e. these orbits are singlets. 

A few orbits in figure ^ do not fit into this simple 
scheme. A closer inspection reveals that these orbits 
are not singlets, but Z-doublets, and indeed they obvi- 
ously occur in pairs. They are generated by symmetry- 
breaking pitchfork bifurcations from singlet orbits. Fig- 
ure [7] presents the orbital parameters for closed orbits 
involved in a bifurcation of this kind. The panels show 
the real and imaginary parts of the starting angles ifi of 
the orbits, the scaled actions S and th e sta bility deter- 
minant M — det(dqf/dpi) (see section [II A ), whose ze- 
ros indicate the occurrence of a bifurcation and that also 
plays an important role in closed-orbit theory . These 
plots should be compared to figure H, which displays the 
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FIG. 7: Orbital parameters close to a pitchfork bifurca- 
tion creating a Z-doublet of closed orbits from a singlet at 
E = —1.4, repetition number fi = 10. tpi is the azimuthal 
starting angle, S the scaled action and M the scaled stabil- 
ity determinant, AS = S — 2tt x 5.898159 was introduced 
for graphical purposes. Thick solid lines: singlet orbit, thin 
solid lines: doublet orbits, dashed line: ghost orbits symmet- 
ric with respect to complex conjugation. Dotted lines indicate 
coordinate axes. 



scenario described by the symmetric cusp catastrophe. 
The qualitative agreement between the catastrophe the- 
ory predictions and the numerical findings is evident. 

A closer look at the asymmetric orbits reveals that they 
have equal initial and final azimuthal angles <p.i — iff, 
i.e. they are not only closed, but also periodic orbits. 
The initial and final angles of these orbits satisfy ip^ — 

—ip^P because they are symmetry partners and tpr 1 ' — 
(2) 

iff because they are periodic. Thus, they must fulfil 

iPi^ = — ¥>j • At the bifurcation, the initial angles of the 
two orbits must coincide, so that a bifurcation can only 
take place when tp+ = or tpi = n, and it actually does 
take place every time one of these conditions is fulfilled. 
This process can be seen in figure ||[ e.g., at S/2ir w 25: 
At F — 0.03, the symmetric orbit has not yet crossed 
the line ifi = n, so that no bifurcation has occurred. At 
F = 0.05, this line has been crossed and two asymmetric 
orbits have been created. 

As the electric field strength is increased, the depen- 
dence of the starting angle on the repetition number 
ceases to be linear. Instead, the curves interpolating 
the functions S(tfi) start to develop humps, so that at 
certain values of S, i.e. at certain repetition numbers, 
more than two possible values of if i exist. This develop- 
ment is illustrated in figure |^. The humps indicate the 
occurrence of tangent bifurcations generating additional 
pairs of singlet orbits. This is the type of bifurcation de- 
scribed by the fold catastrophe ( |l^ ) . Orbital parameters 
for orbits involved in a bifurcation of this kind are shown 
in figure ^. As for the pitchfork bifurcation, a compari- 
son of that figure to the catastrophe theory predictions 
in figure || reveals that the bifurcation is well described 
qualitatively by the fold catastrophe. 

Once additional singlet orbits have been generated in 
a tangent bifurcation, doublet orbits can be generated by 
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FIG. 8: Scaled actions and azimuthal starting angles for pla- 
nar orbits at E = -1.4 and (a) F = 0.2, (b) F = 0.5. Singlets 
are indicated by + symbols, Z-doublets by x symbols. 
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FIG. 10: The bifurcation scenario taking place in the neigh- 
borhood of the tangent bifurcation shown in figure ^. (AS = 
S — 2-7T x 26.2735). Thick solid lines: real singlet orbits, thin 
solid lines: real Z-doublet orbits, dashed lines: symmetric 
ghost orbits, dashed-dotted lines: asymmetric ghost orbits. 
Dotted lines are coordinate axes. 



pitchfork bifurcations in the same way as from the orig- 
inal singlet orbits, i.e. a bifurcation will occur whenever 
a singlet orbit crosses one of the lines cp$ — or ipi = tt. 
This is illustrated in figure [h], which presents the tan- 
gent bifurcation already shown in figure |^ at F m 0.11. 
At F = 0.135, one of the orbits thus generated crosses 
the line (pi = 0, and two doublet orbits are created from 
it. Together, the two bifurcations form what Wang and 
Delos jl8) call the "normal sequence" of bifurcations, 
whereas a pitchfork bifurcation of a singlet orbit gen- 
erated at F = 0, which is not preceded by a tangent 
bifurcation, is called a "truncated series" . These authors 
introduce an integrable model Hamiltonian to explain 
why this kind of sequences can often be observe d for 
planar orbits 
and 
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sheds new light on this question, suggesting 
that normal sequences can actually be expected to occur 
even more generally than surmised by Wang and Delos. 
In particular, although the crossed-fields system is close 
to integrable at the field strengths considered here, inte- 
grability is not needed to make pitchfork bifurcations a 
generic phenomenon. Instead, the presence of a reflection 
symmetry suffices to reduce its codimension from two to 
one. The sequence of a tangent and a pitchfork bifurca- 
tion, represented as a sequence of a fold and a symmetric 
cusp catastrophe, can be regarded as an unfolding of the 
symmetrized version of the butterfly catastrophe [pol |2£ 



FIG. 9: Orbital parameters close to a tangent bifurcation of 
planar orbits at a E = —1.4 and a winding number of fi — 45. 
(AS = S - 2ty x 26.512735.) Solid lines: real orbits, dashed- 
dotted lines: ghost orbits. 
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which is of codimension two, so that its unfolding can be 
expected to occur frequently in codimension one. 

A third bifurcation can be discerned in figure [l^: At 
F w 0.225, a singlet orbit generated at F = and a sin- 
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FIG. 11: Scaled actions and polar starting angles of closed or- 
bits at E — —1.4 and F — (pure magnetic field, + symbols), 
F = 0.05 (x symbols) and F — 0.1 (□ symbols). Due to the 
Z-symmetry, the figure should be extended to be symmetric 
with respect to the line ~&i = 7v/2. 



glet orbit generated in the tangent bifurcation discussed 
above collide and are destroyed. This is an instant of an 
inverse tangent bifurcation, which can be described by 
the fold catastrophe in the same way as the "regular" 
tangent bifurcation. It forms the third building block for 
the bifurcation scenario changing the pattern of planar 
orbits as the electric field strength is increased. 

Besides the three bifurcations described above, in fig- 
ure |l^ three further zeros of the stability determinant M 
occur for certain real orbits, indicating the presence of 
even more bifurcations. These bifurcations involve non- 
planar orbits, i.e. they are pitchfork bifurcations break- 
ing the Z-symmetry all planar orbits possess. They will 
be discussed further in subsequent sections. At the mo- 
ment it suffices to note that in this scenario six individual 
bifurcations take place in a comparatively small interval 
of the electric field strength. This is the first example 
of a phenomenon to be encountered repeatedly: In the 
crossed-ficlds hydrogen atom bifurcations of closed orbits 
abound, exacerbating both the classical and the semiclas- 
sical treatment of the dynamics. 



B. Non-planar orbits 

The transition from the rotationally-symmetric case of 
a pure magnetic field to crossed fields occurs for non- 
planar orbits in much the same way as for planar orbits. 
As soon as a small perpendicular electric field is present, 
a one-parameter family of DKP orbits is destroyed and 
splits into two isolated closed orbits. These orbits start 
in opposite directions with respect to the electric field, 
so that their azimuthal starting angles if i differ by tt, in 
complete analogy with what was shown in figure |. An 
additional complication arises because the polar starting 



angle i?i is no longer bound to the fixed value tt/2, so that 
the two orbits will in general have different di . Figure |ll| 
presents the polar starting angles and the scaled actions 
of the closed orbits for the scaled energy E = — 1.4 in a 
pure magnetic field and for two different perpendicular 
electric field strengths. Only angles $i < tt/2 need to be 
shown because orbits with di > tt/2 can be obtained by 
a Z-reflection. It is obvious from the figure how a family 
of orbits splits in two isolated orbits and the two orbits 
move apart as the electric field strength is increased. This 
process takes place in the same way for both rotator and 
vibrator orbits. 

An exceptional role is played by the DKP orbit parallel 
to the magnetic field. This orbit is isolated even in a pure 
magnetic field. In the presence of a perpendicular electric 
field it is distorted and torn away from the magnetic field 
axis, but it remains isolated rather than splitting into 
two orbits. This process is also apparent from figure [Tl|. 
Notice again that closed orbits in crossed fields do not 
possess repetitions. Any repetition of the parallel DKP 
orbit gives rise to a closed orbit in crossed fields (for 
sufficiently small F), but these orbits are not repetitions 
of each other. They have, in particular, different starting 
angles. 

The symmetries of the closed orbits are worth noting. 
All non-planar orbits described so far are doublets. More 
precisely, the vibrator orbits are T-doublets, i.e. they are 
invariant under the T operation. Their initial and final 
polar angles are small, as the orbits are mainly directed 
along the magnetic field axis. 

For the rotator orbits the situation is more complex. 
Their main component is the motion in the plane perpen- 
dicular to the magnetic field, whereas the z-component 
is comparatively small. They have, therefore, initial and 
final polar angles close to tt/2, so that it is conceivable 
that they can start at an angle di < tt/2 "above" the 
x-y-plane and return at > tt/2 "below" that plane. 
This is in fact the case for the rotators of the first series. 
They turn out to be C-doublets. 

The second series of rotators contains orbits which, in 
the case of a pure magnetic field, are composed of a first- 
series orbit and its Z-reflected counterpart. The orbits 
of the second series therefore have $i = fit and are T- 
doublets. By the same token, orbits of the third series re- 
turn "below" the x-y-plane again and are C-doublets, and 
higher series of rotators alternatingly contain T-doublets 
and C-doublets. 

The distribution of symmetries is illustrated in fig- 
ure |l2](a). It extends the data given in figure [ll] to longer 
orbits and classifies the orbits according to their symme- 
tries. 

So far, only orbits present at arbitrarily low electric 
field strengths have been described. As the electric field 
strength increases, further bifurcations occur. Their gen- 
eral pattern can be identified in figures 11 and [T^(a) . The 
most obvious consequence of the bifurcations is the ap- 
pearance of quartet orbits in each series of both rotator 
and vibrator orbits. They are generated by pitchfork bi- 
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FIG. 12: Scaled actions and polar starting angles of closed 
orbits at E = -1.4 and (a) F = 0.1, (b) F = 0.6. Orbits 
are classified according to their symmetries: T-doublets are 
indicated by + symbols, C-doublets by x symbols, quartets 
by □ symbols. Planar orbits (Z-doublets and singlets) are 
omitted. Notice that the range of actions shown is smaller in 
(b). 
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FIG. 13: Orbital parameters close to a pitchfork bifurcation 
of a first-series rotator and a repetition number of \i = 38. 
The bifurcation creates a quartet of orbits from a C-doublet. 
(AS = S - 2vr x 18.297822.) 




0.08 0.1 0.12 0.14 



0.08 0.1 



0.12 0.14 



FIG. 14: "Normal sequence" of bifurcations for non-planar 
rotator orbits of the second series and a repetition number of 
fx = 54. (AS = S - 2tt x 31.84035.) 



furcations from the adjacent doublet orbits. As figure [l3| 
reveals if it is compared to figure 0, this bifurcation is 
very similar to a pitchfork bifurcation of planar orbits. 
A difference arises because, due to the absence of Z- 
symmetry, the angle $j is not restricted to a fixed value. 
As the C-symmetry to be broken concerns the azimuth 
angles, it still is predominantly the angle ipi that shows 
a square root behavior at the bifurcation and obtains an 
imaginary part when ghost orbits exist. Nevertheless the 
polar angle i?j also acquires a small imaginary part. The 
real part of i?j apparently behaves linear close to the bi- 
furcation, although for electric field strengths above the 
critical value a square root behavior must be present. It 
is too small to be seen in the figure. Even though in quan- 
titative terms the ^-direction is only marginally involved 
in the bifurcation, its presence has the important conse- 
quence that the quartet orbits are no longer constrained 
to be periodic. As the distance from the bifurcation is 
increased, the periodicity condition <pt — ipf is increas- 
ingly, albeit slowly, violated. The same features can be 
found for the bifurcations introducing the quartet orbits 
into the vibrator series. 

The second important type of bifurcations is a tan- 
gent bifurcation introducing new doublet orbits into the 
series. The occurrence of this phenomenon can be no- 
ticed in figure [H], if the numbers of orbits of a given rep- 
etition number are compared for different electric field 
strengths. An example of this bifurcation is given in fig- 
ure [14| The tangent bifurcation involves both angles to 
roughly equal extent. The two doublet orbits thus gener- 
ated are implanted into the regular pattern of their series, 
so that one of them subsequently undergoes a pitchfork 
bifurcation which creates a quartet. This phenomenon is 
entirely analogous to the "normal sequence" of bifurca- 
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FIG. 15: Destruction of T-doublet orbits in a collision with a 
singlet orbit. (AS = S - 2tt x 27.60324.) 
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FIG. 16: Destruction of quartet orbits in a collision with Z- 
doublet orbits. (AS = S - 2tt x 26.569655.) 



tions that was found for planar orbits, except that the 
quartet orbits thus generated are not periodic. 

As the electric field strength increases, the rotator or- 
bits of a given series are torn apart and span an ever 
wider interval of Those orbits moving towards higher 
values of #j eventually hit the plane di = tc/2, where they 
collide with their Z-refiected partners and are destroyed. 
One might suspect the destruction of the two orbits to 
occur in a tang ent bifurcation, but from the discussion 
of section 
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it is clear that a tangent bifurcation can 
only create or destroy orbits having different actions, so 
that it can never involve two orbits related by a sym- 
metry transformation. Thus, the bifurcation must be of 
pitchfork type, and it must involve a Z-symmetric planar 
orbit. Depending on whether the non-planar orbits col- 
liding with the plane are doublets or quartets, the planar 
orbit must be a singlet or a Z-doublet, respectively. If 
the destruction scenario is regarded in the direction of 
decreasing field strengths, it appears as the creation of 
orbits with broken Z-symmetry from an orbit possessing 
this symmetry. It is therefore the Z-breaking analogue of 
the T- and C-symmetry breaking bifurcations described 
above. As this type of bifurcation involves a planar orbit, 
it must give rise to a zero in the stability determinant M 
of t he p lanar orbit. In fact, the examples given in fig- 
ures [Tq and [l^ for both the destruction of a doublet and 
a quartet are two of the three bifurcations whose presence 
was inferred from figure [l^ the discussion of the planar 
orbits. 
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FIG. 17: Simple bifurcation scenario for vibrator orbits of 
repetition number fx = 41. (AS = S - 2tt x 24.50221.) 



The scenario just described is not restricted to rotator 
orbits. As can be seen in figure ll], the short vibrator or- 
bits can, even at low electric field strengths, reach rather 
high values of At F — 0.15550 the first of them col- 
lides, at "&i = 7r/2, with its Z-refiected counterpart and 
is annihilated. This is a pitchfork bifurcation in which 
one of the planar orbits with repetition number fj, = 1 
takes part. Similarly, longer vibrators are destroyed in 
collisions with planar rotators of the appropriate repeti- 
tion numbers. This example demonstrates that the dis- 
tinction between vibrators and rotators, which was bor- 
rowed from the case of vanishing electric field, does not 
apply, strictly speaking, if an electric field is present. Al- 
though it is generally useful for rather high electric field 
strengths, it can fail in some instances. This is clearly 
the case when a bifurcation involves both vibrator and 
rotator orbits. 

A collision with the plane perpendicular to the mag- 
netic field occurs only for vibrators of low repetition num- 
bers, and only for vibrators that descend from the orbit 
parallel to the magnetic field. For longer orbits, the usual 
scenario is different. At low electric field strength there is, 
for sufficiently high repetition numbers, one orbit stem- 
ming from the orbit parallel to the magnetic field and 
one or several pairs of orbits created from non-parallel 
vibrators. It can be seen in figure [ll], however, that 
for certain repetition numbers two of these orbits can 
be missing. This happens when the descendant of the 
parallel orbit and one of the other vibrators annihilate 
in a tangent bifurcation. A simple example of how this 
can come about is provided by the orbits with the repe- 
tition number /i = 41. Their bifurcations are illustrated 
in figure [l^. Two of the orbits obviously bifurcate from 
a common family at F = 0, whereas the orbit proceed- 
ing from the parallel orbit is isolated there and starts at 
di = 0. It then merges with one of the other orbits in a 
tangent bifurcation to form a pair of ghost orbits. 

This bifurcation is as simple as one could expect. For 
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FIG. 18: Complicated bifurcation scenario for vibrator orbits 
of repetition number /j = 42. (AS = S - 2tt x 25.09941.) 



the neighboring vibration number \i — 42 the scenario is 
more complicated. It is illustrated in figure [l|. In this 
case, one of the orbits generated in the rotational symme- 
try breaking at F = 0, which is a T-doublet, undergoes 
a pitchfork bifurcation and gives birth to a quartet of 
orbits before it annihilates with the descendant of the 
parallel orbit. The quartet then collides with the third, 
leftover T-doublet and is destroyed in a second pitchfork 
bifurcation. 

Corresponding to the three elementary bifurcations, 
there are three ghost orbits involved in the scenario. For 
one of them, the starting angles $j and ip* show a pecu- 
liar behavior at the electric field strength Fq = 0.08750: 
Whereas i?j exhibits a square root behavior, changing 
from nearly real to nearly imaginary values, the real part 
of tfi changes discontinuously by tt/2, and the imaginary 
part of <pi seems to diverge. Neither the action nor the 
stability determinant of the orbit, on the contrary, show 
any kind of special behavior. In particular, M is non- 
zero, so that there cannot be a bifurcation of the ghost 
orbit. 

The Cartesian components of the unit vector s in 
the starting direction are given in figure |l9|. For all of 
them either the real or the imaginary parts are small, 
so that their numerical calculation is hard. Neverthe- 
less, to within the numerical accuracy all components 
are smooth at Fq, although the angles i?j and tpi used 
to calculate them are not. Thus, the singularity must 
be due to the transformation from Cartesian to angular 
coordinates. In the real case it is obvious that the <p) 
coordinate chart is singular at i? = 0. To elucidate the 
details in the case of ghost orbits, assume a model sit- 
uation where s z = cos is exactly real and s z = 1 at 
F = Fq. For ghost orbits, s z is not bound to be smaller 
than 1, so that generically, to first order in e = F — Fq, 
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FIG. 19: Cartesian components of the unit vector s specifying 
the starting direction of the ghost orbit (s x — sin i?i cos ipi , 
s y = sin$i sin ipi, s z — cosi?j). Vertical dotted lines mark the 
field strength Fq — 0.08750 where the singularity of Im cpi is 
encountered. 

cosf?i — 1 cx e. Therefore, cx y/E shows a square root 
behavior and changes from purely real to purely imagi- 
nary values. At the same time, sini?j cx \fe has a zero, 
so that for s x = sin i?j cos tpi and s y = sin $j sin tpi to 
assume finite values, sin ipi and cos tpi must diverge as 
£ -i/2 This is only possible if the imaginary part of tp-i 
is large. More precisely, if Im ipi > is assumed to be 
large, sincpi and cost/?i are proportional to e~ l<Pi , whence 

Vi = ^\ne + O(e ) (22) 

achieves the desired divergence of sin tpi and cos tfi . Now 
Re lne = In |e| diverges at e = 0, whereas Im lne changes 
discontinuously from to ±7r, depending on what branch 
of the logarithm is chosen. This behavior results in the 
observed divergence of Im c/?j and a discontinuous jump 
in Re (/Si of size tt/2. 

In the actual scenario s z will not be exactly equal to 1 
at Fq because this is a situation of real codimension two. 
However, if Im s z is small, the singular behavior described 
above will be closely approximated. Indeed, a closer look 
at i?j (see figure ^0|) reveals that it is actually smooth, but 
close to F — Fq it changes extremely rapidly from almost 
real to predominantly imaginary values. Similarly, the 
real part of tpi is smooth, although it changes over an 
even smaller range of F. From the numerical data it 
cannot be determined if Im tpi is also smooth or actually 
diverges. From the above discussion it is clear that it 
must be smooth, because the coordinate singularity at 
•di = is not actually encountered. 

It should be noted that the singularity described here 
can occur for ghost orbits only. In the real case, as the 
pole $i — on the real unit sphere (which still has codi- 
mension two) is approached, both s x and s y must vanish 
instead of assuming finite values, so that no divergences 
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FIG. 20: The starting angle i?i of the ghost orbit close to 
the coordinate singularity. Solid line: real part, dashed line: 
imaginary part. Note the logarithmic scale. 



of any kind are required. 



VI. THE CLASSIFICATION OF CLOSED 
ORBITS 

The fundamental classification scheme used in the 
above description of closed orbit bifurcations is the dis- 
tinction between rotators and vibrators. This distinc- 
tion was adopted from the case of vanishing electric field 
strength, so it can be expected to be applicable if the 
electric field is not too strong. In fact, from figure |l^(a) 
it is obvious that for E — —1.4 and F = 0.1 orbits can 
uniquely be classified as rotators or vibrators and can be 
assigned both a series number and a repetition number. 
However, for long orbits neighboring series of rotators 
already start to overlap in the plot, and if the electric 
field strength is increased to F = 0.6, all orbits get com- 
pletely mixed up, producing the somewhat messy picture 
shown in figure [l2](b) . The relative strengths of the elec- 
tric and magnetic fields can conveniently be compared if, 
instead of scaling the magnetic field strength to 1 as de- 
scribed in section ||, the scaling properties of the Hamil- 
tonian are used to scale the energy to E = — 1. For the 
present parameter values, the scaled field strengths are 
then B = 0.60 and F = 0.30, so that the magnetic field 
is still stronger than the electric field, although the latter 
is no longer negligibly small. Thus, the pattern of or- 
bits should still be dominated by the magnetic field, but 
figure |l2](b) suggests that there is no way to achieve a 
classification. 

A more suitable starting point for a classification is 
provided by the complete trajectories. Of course, since 
the classification must gradually break down for suffi- 
ciently strong electric fields, it can only be based on 
heuristic criteria. The criteria we are going to propose 
are largely based on the behavior of the z-coordinate of 
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FIG. 21: Rotor orbits of the (a) first, (b) second and (c) 
third series: scaled coordinate z as a function of the scaled 
pseudotime f for E = —1.4 and F = 0.2. 



the motion as a function of time. To illustrate them, this 
function is plotted for rotators of roughly equal length 
from different series in figure [2l]. Figure shows the 
analogous data for vibrators. 

First of all, vibrators are connected to an orbit along 
the z-axis in the pure magnetic field case. In this limit, 
the motion takes place either "above" the x-y-plane, i.e. 
in the half-space z > 0, or "below" the plane. Rota- 
tors, on the contrary, stem from the elementary orbit in 
the plane. Their motion takes place both above or be- 
low the plane. Rotators can therefore be distinguished 
from vibrators if the maximum and minimum values of 
the coordinate z are compared: For a rotator, they must 
have roughly equal absolute values, whereas for a vi- 
brator "above" the plane, the minimum value is much 
smaller in magnitude than the maximum value. 

For the vibrators of a given length shown in figure |2^, 
this criterion gets better the higher the series of the vibra- 
tor is chosen. For the vibrator of the first series, which 
is closest to the domain of rotators, the excursion into 
the lower half space is of the same order of magnitude 
as that into the upper half space. As the electric field 
strength increases further, the vibrator orbit will even- 
tually become indistinguishable, by the present criteria, 
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FIG. 22: Vibrator orbits of the (a) first, (b) second and (c) 
third series: scaled coordinate z as a function of the scaled 
pseudotime f for E = —1.4 and F = 0.2. 



from a rotator of the second series. 

It has already been noted in the discussion of orbital 
symmetries that a rotator of the first series that starts 
from the nucleus into the upper half space returns to it 
from the lower half space, whereas a rotator of the second 
series returns from the upper half space. This alternation 
between motion in the upper and lower half spaces is 
obvious from figure ^l|. It can be used to determine the 
series of a rotator. If the value of z{t) in a maximum is 
compared to its value in the subsequent minimum, the 
transitions between motion in the upper and lower half 
spaces can easily be monitored. 

Assigning a series number to a vibrator is considerably 
more difficult. It relies on the beat structure present in 
z(t) for a vibrator. Subsequent maxima of this func- 
tion have varying heights. We consider the maxima that 
have minimal height compared to neighboring maxima, 
i.e. the minima of the beats apparent in figure ^2|. The 
height of these beat minima slowly oscillates, and the 
number of their oscillations turns out to give the series 
number of the vibrator. It can be found by counting 
the number of minima in the oscillation. Notice that the 
method uses extrema of the third order: the minima in 
the minima in the height of maxima. It therefore requires 
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FIG. 23: Rotators of the first series for E — —1.4 and (a) 
F = 0.2, (b) F = 0.4, and (c) F = 0.6. C-doublets are 
indicated by x symbols, quartets by * symbols. 



vibrators of sufficiently high repetition numbers, so that 
many maxima of z(t) exist. The method will fail when 
applied to vibrators of very small repetition numbers. 
However, these vibrators exist for fairly high scaled en- 
ergies only, where even in the absence of an electric field 
the dynamics is mixed or chaotic. In these regions, the 
classification suggested here might not be meaningful at 
all. 

Finally, orbits can be assigned a repetition number. 
For vibrators, this can be done by simply counting the 
number of maxima in z(t). For rotators, maxima in p(r) 
must be counted, where p — {x 2 + y 2 ) 1 ^ 2 is the distance 
from the magnetic field axis. In this case, an additional 
difficulty arises because p cannot assume negative values, 
so that between two maxima corresponding to repetitions 
of the elementary orbits, additional maxima arise from 
the "swing- by" around the nucleus. These maxima tend 
to be extremely small and narrow close to the beginning 
and the end of an orbit, so that they are hard to detect, 
but they can reach considerable heights in the middle 
of an orbit. The present form of the classification algo- 
rithm, which is certainly open to improvements, counts a 
maximum in p(r) if its height exceeds a certain fraction 
of the height of the previous maximum. In this form, the 
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assignment of rotator repetition numbers turns out to be 
the least robust part of the classification procedure. 

The criteria described above readily lend themselves 
to a numerical implementation, so that the classification 
of orbits can be achieved automatically. As an exam- 
ple, the rotators of the first series are shown in figure |23| 
for three different electric field strengths. Although the 
neat "wiggly-line" structure characterizing the series in 
figure [ll] quickly breaks down for larger electric field 
strengths, the distinction between different series per- 
sists. Figure ^3|(c) should be compared to figure |l^(b). 
ft might appear surprising that the messy-looking set of 
orbits still permits a classification, but with the help of 
the criteria just described an ordered pattern of closed 
orbits can still be discerned. In this sense, the classi- 
fication scheme derived from the DKP turns out to be 
remarkably robust. 

VII. SUMMARY 

In this work, a systematic study of the closed classi- 
cal orbits of the hydrogen atom in crossed electric and 



magnetic fields has been carried out. As an important 
step towards a complete understanding of the compli- 
cated pattern of closed orbits, a bifurcation theory of 
closed orbits was developed and the generic bifurcations 
of codimension one were identified. 

A variety of bifurcation scenarios observed in the 
crossed-fields system was described. They demonstrate 
that, even though only two types of elementary bifurca- 
tions exist, they combine into a wealth of complicated 
bifurcation scenarios. The abundance of bifurcations ex- 
acerbates both a complete classical description of the 
crossed-fields hydrogen atom and its semiclassical treat- 
ment 

Based on the classification of closed orbits in the hydro- 
gen atom in a magnetic field, heuristic criteria have been 
proposed which allow a systematization of closed orbits 
for moderately high electric field strengths. Although 
the present analysis cannot yet claim to have achieved 
a complete classification of closed orbits in the crossed- 
fields hydrogen atom, it does give a detailed impression 
of how orbits bifurcate as the electric field strength in- 
creases. It thus introduces a high degree of order into the 
complex set of closed orbits. 
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